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SUMMARY 


The  shock  Hugoniot  and  adiabatic  release  curves  of  sedimentary  rocks 
display  a  range  of  different  behaviors  that  affect  the  decay  of  shock  waves 
propagating  away  form  a  confined  source.  Hysteresis  in  the  shock-release  paths  of 
geologic  materials  results  in  irreversible  energy  deposition,  thus  depriving  the 
shock  wave  of  energy  required  to  drive  its  propagation  and  resulting  in  a  more  rapid 
wave  amplitude  decay  and  less  efficient  coupling  to  far  field  seismic  waves  than 
would  be  expected  from  geometrical  effects  alone.  Thus,  data  constraining  this 
behavior  in  rocks  are  essential  for  determining  the  cavity  volumes  for  various 
degrees  of  decoupling  for  different  lithologies. 

We  have  studied  the  shock  compression  and  adiabatic  release  behavior  of 
Coconino  Sandstone  with  an  initial  density  p^o  =  2.313  Mg/m^  and  Solenhofen 
Limestone  with  initial  density  Poo  =  2.613  Mg/m^,  a  Hmestone  from  Bedford, 
Indiana,  with  initial  density  Pog  =  2.418  Mg/m^  and  a  calcareous  shale  from  the 
Chambless  formation  with  initial  density  poo=  2.576  Mg/m^.  We  used  a  series  of 
shock  Hugoniot  and  release  adiabat  measurements  on  sandstone,  limestone,  and 
shale  to  develop  models  for  the  shock  and  release  behavior  of  these  rocks.  The 
sandstone  shock  Hugoniot  curve  can  be  calculated  from  the  equations  of  state  of 
quartz,  stishovite,  and  a  high  coordination  number  liquid.  A  new  fit  of  available 
data  for  the  dense  liquid  phase  gives  poo  =  4.150  Mg/m^,  Kso  =  313.5  GPa,  K'  = 
1.8484,  y  =  1.4(V/Vo)°•^  and  En  =  2.4  MJ/kg.  In  the  case  of  Coconino 
Sandstone,  the  presence  of  a  few  percent  of  non-quartz  constituents  requires  that 
the  STP  densities  of  the  relevant  phases  be  reduced  by  2.4%  to  fit  the  data.  The 
mixed-phase  regions  of  the  Hugoniot  can  be  modeled  as  a  simple  mixture  of  the 
two  relevant  end  member  phases.  Formation  of  the  higher  pressure  phase  appears 
to  be  kinetically  inhibited,  with  the  fraction  f  of  the  higher  pressure  phase  being 
described  by  f  =  Ae"^*^®’^ .  For  formation  of  stishovite  from  quartz,  A  =  7.286  and 
Ea(kJ  /  mol)  =  89.36  -  71.97(po/Poo )  where  poo  and  po  are  the  actual  and 

theoretical  nonporous  densities  of  the  initial  rock  sample.  For  formation  of  the 
dense  liquid  from  stishovite,  A  =  6.245  x  10^  and  E^  (kJ/mol)  =  404  +  0.312poo. 

No  formation  of  the  higher  pressure  phase  occurs  until  the  phase  boundary  is 
reached.  As  crude  approximations  we  describe  the  quartz-stishovite  phase 
boundary  by  P  =  12  GPa  and  the  stishovite-dense  liquid  phase  boundary  by 
T  =  4348  -  12.121P  for  T  in  K  and  P  in  GPa. 


Release  adiabats  from  the  single  phase  regions  of  the  sandstone  Hugoniot 
are  modeled  by  the  isentropes  of  the  relevant  phases.  At  ~5  GPa,  stishovite  begins 
to  form  a  diaplectic  glass.  We  model  the  release  adiabat  below  this  pressure  as  a 
straight  line  in  P-V  space,  from  the  stishovite  release  isentrope  at  5  GPa  to  the  STP 
diaplectic  glass  density  of  2.27  Mg/m^.  The  dense  liquid  undergoes  a  similar 
change  beginning  at  ~5  GPa  with  formation  of  normal  silica  glass  with  an  STP 
density  of  2.204  Mg/m^.  The  quartz  release  adiabat  follows  the  quartz  isentrope  to 
P=0.  As  with  the  Hugoniot,  we  model  the  release  from  the  mixed-phase  regions  of 
the  Hugoniot  as  a  mixture  of  the  relevant  phases,  with  the  proportions  frozen  at  the 
Hugoniot  values. 

The  limestone  Hugoniot  curve  is  complicated  by  several  metastable  phases 
and  multiple-wave  shock  structure  at  low  pressures.  Single  phase  portions  of  the 
Hugoniots  of  Solenhofen  and  Bedford  limestones  are  modeled  as  calcite  I  at  low 
pressures  and  a  dense  liquid  at  high  pressures.  We  have  fit  an  equation  of  state  for 
the  dense  liquid  phase,  with  po  =  3.377  Mg/m^  Kso  =139  GPa,  K’  =  3.4, 7  = 
2.5(VA^o)i  and  E,,  =  2.142  MJ/kg.  At  intermediate  pressures,  dense  limestones 
such  as  Solenhofen  seem  to  form  calcite  ni,  while  the  more  porous  Bedford 
limestone  apparently  forms  aragonite.  The  mixed  phase  regions  of  the  Hugoniot 
are  modeled  in  a  manner  analogous  to  that  used  for  the  sandstone  Hugoniot.  The 
parameters  describing  the  fraction  f  of  the  higher  pressure  phase  in  the  mixed  phase 
regions  are  A  =  5.848  and  E,  =  62.71 3-20.65  Ipoo  for  formation  of  calcite  IB  or 
aragonite  from  calcite  I  and  A  =  2.919  and  Ea  =  382.453-1 30. 448poo  for  formation 
of  the  dense  liquid  from  either  aragonite  or  calcite  IB,  for  Eg  in  kJ/mol  and  poo  in 
Mg/m\  The  phase  boundaries  in  our  model  are  approximated  by  Tj.^r  =  164.7P- 
3.6,  Ti.n,  =  78.00-1.3,  T^^uq  =  745+48P,  and  Ti„.uq  =  401+25.8P. 

As  with  sandstone,  we  model  the  release  adiabats  with  the  isentropes  of  the 
relevant  phases,  mixing  them  at  the  appropriate  ratios  for  release  from  the  mixed- 
phase  portions  of  the  Hugoniot.  There  is  no  evidence  of  formation  of  a  glass  at 
very  low  pressures,  but  release  of  phases  other  than  calcite  I  is  apparently 
accompanied  by  decomposition  with  evolution  of  CO2  gas  at  pressures  below  1 
GPa.  We  have  ignored  this  very  low  pressure  effect  in  our  calculations. 

Equation  of  state  parameters  for  shale  are  calculated  from  the  mineral 
composition,  with  the  bulk  modulus  and  its  pressure  derivative  constrained  by  the 
shock  wave  data.  The  shale,  which  is  composed  of  roughly  equal  amounts  of 
quartz,  calcite,  and  phyllosilicates,  has  a  low  pressure  phase  characterized  by  po= 
2.615  Mg/m3  K,„  =  30.9  GPa,  K'  =  8.15,  and  7  =  0.73  (v/vo)>  ^2,  and  a  high- 
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pressure  phase  with  po  =  3.659  Mgkm^,  K^o  =  265  GPa,  K'  =  3.14,  and  y=  1.2 
(VA^o)*'^®-  The  release  adiabats  calculated  from  these  equations  of  state  generally 
agree  with  the  data,  but  do  not  take  into  account  large-volume  phase  changes  that 
occur  at  low  pressure. 

These  results  allow  us  to  estimate  the  fraction  f  of  shock  internal  energy 
irreversibly  deposited  in  these  rocks  and  thus  no  longer  available  to  drive  the  shock 
wave.  We  find  that  f  is  a  factor  of  -2  smaller  in  limestones  than  in  sandstone,  so 
that  decoupling  of  the  source  from  far-field  seismic  waves  is  more  efficient  in 
sandstones.  The  value  of  f  may  be  reduced  considerably  in  the  presence  of  an 
elastic  precursor  wave  and  of  intermediate  waves.  Experimental  evidence  indicates 
that  shock  wave  decay  still  occurs  even  in  the  supposedly  purely  elastic  shock 
regime,  suggesting  that  effects  such  as  fracture  formation  are  very  important  to  the 
overall  picture  of  explosive  energy  coupling. 


L  INTRODUCTION 


The  shock  Hugoniot  and  adiabatic  release  curves  of  sedimentary  rocks 
display  a  range  of  different  behaviors  that  affect  the  decay  of  shock  waves 
propagating  away  form  a  confined  source.  Many  of  the  minerals  making  up 
sandstones,  shales,  and  limestones  undergo  phase  transitions  when  shocked  to 
pressures  of  interest  to  studies  of  coupling  of  energy  from  explosive  sources  into 
far  field  seismic  waves.  Both  sandstones  and  limestones  have  been  observed  to 
exhibit  elastic  precursors  and  multiple  wave  behavior  at  shock  velocities  up  to  3.7 
km/s  and  5.7  km/s,  respectively  [Ahrens  and  Gregson,  1964].  Hysteresis  in  the 
shock-release  paths  of  materials  results  in  irreversible  energy  deposition,  thus 
depriving  the  shock  wave  of  energy  required  to  drive  its  propagation  and  resulting 
in  a  more  rapid  wave  amplitude  decay  and  less  efficient  coupling  to  far  field  seismic 
waves  than  would  be  expected  from  geometrical  effects  alone.  Thus,  data 
constraining  this  behavior  in  rocks  are  essential  for  determining  the  cavity  volumes 
for  various  degrees  of  decoupling  for  different  lithologies. 

We  have  studied  the  shock  compression  and  adiabatic  release  behavior  of 
Coconino  Sandstone  with  a  mean  initial  density  of  =  2.313  Mg/m^  and 
Solenhofen  Limestone  with  a  mean  initial  density  of  Pqq  =  2.613  Mg/m^,  a 
limestone  from  Bedford,  Indiana,  with  a  mean  initial  density  of  Poo  =  2.418  Mg/m^, 
and  a  2.62  Mg/m^  calcareous  shale  (Chambless  formation).  The  sandstone  is 
primarily  quartz,  with  ~3%  feldspar.  The  limestones  are  essentially  pure  calcite 
with  very  minor  amounts  of  other  constituents.  The  shale  is  nearly  equal  fractions 
of  calcite,  clays  and  micas,  and  quartz. 
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IL  SHOCK  WAVE  EQUATION  OF  STATE  OF  SANDSTONE, 
LIMESTONE,  AND  SHALE 


A.  Experimental  Techniques. 

The  experimental  techniques  used  for  this  study  are  discussed  in  detail 
elsewhere  [Ahrens,  1987].  A  sample  of  the  rock  being  studied  is  placed  on  a  metal 
driver  plate  (Figure  1),  with  standard  buffers  and/or  inclined  mirrors  placed  on  the 
free  (rear)  surface  of  the  sample.  Shock  arrival  mirrors  are  placed  against  the  rear 
surfaces  of  the  driver  plate,  sample,  and  buffers.  During  the  experiment,  the  driver 
plate  is  impacted  by  a  gun-launched  projectile  containing  a  metal  impactor  (flyer 
plate).  The  impact  generates  a  shock  wave  which  propagates  through  the  driver 
plate,  sample,  and  buffer  materials.  Arrival  of  the  shock  wave  at  the  free  surface  of 
each  component  of  the  target  causes  disruption  of  the  specular  reflection  from  the 
shock  arrival  mirrors  as  shown  in  Figure  1  b.  One  spatial  dimension  of  the  target, 
including  all  shock  arrival  mirrors,  is  imaged  as  a  function  of  time  by  a  streak 
camera,  allowing  the  shock  wave  transit  times  and  thus  the  shock  wave  velocities 
through  target  components  of  known  thickness  to  be  determined.  The  internal 
energy  E,  pressure  P,  density  p,  and  particle  velocity  Up  of  the  shocked  sample  and 
buffer  are  calculated  from  the  measured  projectile  velocity,  sample  initial  density, 
and  shock  wave  velocity  using  the  Rankine-Hugoniot  shock  equations  for 
conservation  of  mass,  momentum,  and  energy  and  by  the  requirement  that  pressure 
and  particle  velocity  must  be  continuous  across  material  interfaces. 

The  buffers  and  inclined  mirrors  are  used  to  measure  points  on  the  release 
adiabats  of  the  samples  by  the  requirement  that  the  P-Up  state  in  the  released  state  of 
the  sample  be  identical  with  that  in  the  shocked  buffer,  or  the  free  surface  in  the 
case  of  the  inclined  mirror.  Table  1  give  the  properties  of  the  materials  used  for 
flyer  and  driver  plates  and  buffers  in  this  study. 

The  shock  Hugoniot  of  a  material,  which  is  the  locus  of  final  states  achieved 
by  shock  compression,  is  conveniently  described  as  a  polynomial  representation  of 
the  shock  wave  velocity,  in  terms  of  the  particle  velocity,  Upi 

U,=Co+sUp-hs’Up'  (1) 
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Figure  1.  Diagram  of  a  typical  shock  Hugoniot  experiment  with  buffers  for  release 
adiabat  measurements,  (a)  Projectile  in  flight  prior  to  impact  (b)  Static  (left)  and 
dynamic  image  recorded  by  streak  camera. 
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Table  1.  Shock  Hugoniot  parameters  for  standard  materials  used  in  this  study. 


Material 

po  (Mg/m3) 

Co  (km/s) 

s 

s'  (s/km) 

Source 

All  100 

2.707 

5.386 

1.339 

1 

A12024 

2.784 

5.330 

1.34 

2 

W 

19.224 

4.029 

1.237 

1 

Ta 

16.650 

3.293 

1.308 

2 

Lexan 

1.193 

2.449 

1.498 

3 

4.817 

0.589 

4 

1.907 

1.431 

5 

Polystyrene 

0.0474 

0.362 

0.9165 

0.04982 

6 

(1)  Marsh  [1980] 

(2)  Mitchell  and  Nellis  [1981] 

(3)  Fit  to  data  from  Marsh  [1980]  for  Up  <  2.605km  /  s. 

(4)  Fit  to  data  from  Marsh  [1980]  for  2.605km  /  s  <  Up  <  3.456km  /  s. 

(5)  Fit  to  data  from  Marsh  [1980]  for  Up  >  3.456km  /  s. 

(6)  Calculated  from  Poo  using  expression  from  Anderson  and  Ahrens  [1993] 
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Usually,  the  polynomial  is  truncated  at  first  or  second  order.  The  pressure  P, 
density  p,  and  internal  energy  E  in  the  shocked  state  are  given  by 


~  PoflUsUp 


(2) 


Ph  “  Poo 


U.-Up 


(3) 


(4) 


Here,  the  subscripts  0  and  H  refer  to  the  unshocked  and  shocked  (Hugoniot)  states, 
respectively,  and  the  subscript  00  on  p  refers  to  the  initial  density,  including  any 
porosity  in  the  sample.  We  arbitrarily  set  E  =  0  as  the  initial  condition. 

B.  Sandstone 

1.  Experimental  Data 

Table  2  presents  the  experimental  data  for  the  Coconino  Sandstone  used  in 
this  study,  which  had  a  mean  initial  bulk  density  of  2.313  Mg/m^.  The  densities  of 
the  release  states  are  upper  limits,  obtained  by  approximating  the  release  paths  as 
straight  line  segments  in  the  P-V  plane  [Lyzenga  and  Ahrens,  1978]. 

Figure  2  shows  the  present  data,  along  with  an  empirical  fit  to  the  Uj-Up 
projection  of  the  Hugoniot.  The  U,-Up  projection  of  the  Hugoniot  curve  can  be 
divided  into  two  straight-line  segments  with  weighted  fits  of  Cq  =  3.210±.026  km/s 
and  s  =  0.978+.01 1  for  Up  <  2.932  km/s  and  =  0.973±.108  km/s  and  s  = 
1.741±.036  for  Up  >  2.939  km/s.  Also  shown  is  the  P-V  projection  of  the  curve 
described  by  the  empirical  U^-Up  fit,  along  with  the  present  data  and  data  from  Van 
Thiel  [1977]  for  comparison. 

Sandstones  generally  show  an  elastic  precursor  for  very  low  shock 
pressures.  Based  on  the  data  from  Ahrens  and  Gregson  [1964],  the  velocity  Ug,  of 
the  elastic  precursor  is  only  very  weakly  dependent  on  the  particle  velocity,  with 
even  the  sign  of  the  variation  being  uncertain.  As  a  first  order  approximation,  the 
elastic  precursor  velocity  depends  only  on  the  initial  density  of  the  rock  .  Using  the 
data  from  Ahrens  and  Gregson  [1964],  which  cover  polycrystalline  quartz  rocks 
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Table  2.  Shock  wave  and  release  adiabat  data  for  Coconino  Sandstone. 
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Figure  2.  Shock  Hugoniot  and  release  adiabat  data  and  Hugoniot  based  on 
empirical  U.-Up  fits  for  Coconino  Sandstone,  (a)  U,-Up  plane,  (b)  P-V  plane. 


7 


TJA95083SFD 


with  Poo  in  the  range  from  1.961  Mg/m^  to  2.650  Mg/m^,  we  describe  the  elastic 
precursor  velocity  Uei  by 

U„=4.53poo -6.122  (5) 

where  Ug,  is  given  in  km/s  and  Poo  is  in  Mg/m^. 

2.  Equation  of  State. 

Although  the  experimental  U^-Up  relation  gives  us  a  good  empirical 
description  of  the  Hugoniot  of  the  particular  rock  we  studied,  we  wish  to  develop  a 
model  for  shock  and  release  behavior  that  is  generally  applicable  to  quartz  rocks  in 
general.  Such  a  model  could  be  easily  modified  to  account  for  minor  non-quartz 
constituents. 

Since  the  predominant  constituent  is  quartz,  the  phases  appearing  on  the 
quartz  shock  Hugoniot  curve  are  expected  to  also  appear  on  the  sandstone  Hugoniot 
curve.  In  particular,  these  are  quartz,  stishovite,  and  dense  Si02  liquid.  Above  the 
Hugoniot  elastic  limit  (HEL),  pressure  at  a  given  specific  volume  =  1/pH  on  the 
Hugoniot  can  be  calculated  by 

P  —  ^tr  ~  Ps^h/Th  (f.) 

”"(Voo-Vh)/2-Vh/Th 

where  is  the  shock  state  (Hugoniot)  pressure,  Pj  is  the  pressure  of  the  material 
on  the  principal  isentrope  at  V  =  Vh,  Vh  is  the  shock  state  specific  volume,  Vqo  is 
the  initial  specific  volume  of  the  material,  E^,  is  the  phase  transition  energy  at  STP, 
Eg  is  the  internal  energy  gained  in  compression  of  the  material  along  the  principal 
isentrope  to  specific  volume  Vh,  and  Yh  =  Gruneisen  parameter  of  the  material  at 
volume  Vh.  This  expression  equates  the  energy  at  pressure  and  specific  volume 
Vh  obtained  by  following  the  thermodynamic  path  of  the  shock  compression  (a 
straight  line  in  P-V  space)  with  that  obtained  via  phase  change  (if  needed)  at  room 
temperature  and  1  bar  pressure,  isentropic  compression  to  volume  Vh,  and 
isochoric  heating  to  pressure  Py.  We  assume  that  the  principal  isentrope  is 
described  by  the  third  order  Eulerian  finite  strain  Birch-Murnaghan  equation: 
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(8) 


fV  9 

’x" 

x2 

1  . 

- +  --t 

4 

2 

4 

,6  4  2  61 

where 

x  =  {VJV)“’  (9) 

and 

^  =  1(4-K')  (10) 

4 

where  Kjo  Is  the  value  of  the  isentropic  bulk  modulus  Kg  at  STP  and  K’  in  the 
isentropic  pressure  derivative  of  Kg  at  STP.  We  also  assume  that  the  Griineisen 
parameter  y  obeys 

Y  =  Y«(V/V.)'  (11) 

For  the  properties  of  quartz,  we  use  the  equation  of  state  (EOS)  parameters 
obtained  by  averaging  the  values  given  by  Sumino  and  Anderson  [1984].  For 
simplicity,  we  assume  that  Cv=31c/atom.  For  stishovite,  we  use  the  EOS 
parameters  tabulated  by  Lyzenga  et  al.  [1983],  but  with  a  slightly  different  value  of 
Kgo,  which  seems  to  better  describe  the  shock  wave  compression  data  for  quartz 
given  by  Wackerle  [1962]  and  presented  graphically  by  Chhabildas  [1986]  and 
Lyzenga  et  al.  [1983].  The  shock  temperature  data  of  Lyzenga  et  al.  [1983]  suggest 
that  Si02  is  liquid  above  1 17  GPa  and  70  GPa  on  the  quartz  and  fused  quartz  shock 

Hugoniots,  respectively.  Since  our  data  are  for  samples  which  are  only  slightly 
more  dense  than  fused  quartz  and  the  Coconino  Sandstone  data  tabulated  by  Van 
Thiel  [1977]  are  for  densities  lower  than  that  of  fused  quartz,  we  can  expect  our 
highest  pressure  data  to  be  for  shock  states  in  the  liquid  phase.  The  densities  of  the 
shock  states  for  quartz  and  fused  quartz  in  the  relevant  pressure  ranges  are 
incompatible  with  the  known  properties  of  liquid  Si02  at  low  pressure  and  strongly 

suggest  that  the  liquid  phase  on  the  Hugoniot  is  characterized  by  in  6-fold 
coordination  with  O^',  as  it  is  in  stishovite.  This  is  a  reasonable  conclusion,  since 
stishovite  is  the  solid  which  would  undergo  equilibrium  melting  to  give  this  liquid. 
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Lyzenga  et  al.  [1983]  give  a  set  of  EOS  parameters  for  the  dense  liquid  phase,  but 
we  cannot  reconcile  the  resulting  shock  compression  curves  with  all  the  existing 
Hugoniot  data  for  quartz  and  fused  quartz.  We  have  refit  the  values  of  Kjo  and  K' 
to  data  from  Lyzenga  et  al.  [1983]  for  quartz,  Trunin  et  al.  [1971a]  for  porous 
quartz,  Trunin  et  al.  [1971b]  for  quartz  and  quartzite  and,  Lyzenga  et  al.  [1983]  and 
Marsh  [1980]  for  fused  quartz.  We  retain  the  empirical  expression  for  Cy  given  by 
Lyzenga  et  al,  [1983].  The  resulting  fit  is  shown  in  Figure  3.  The  EOS  parameters 
used  for  the  three  phases  of  interest  to  this  study  are  listed  in  Table  3,  along  with 
slight  variations  of  these  values  which  give  a  better  fit  to  the  data  for  Coconino 
Sandstone.  The  primary  change  made  to  bring  the  model  in  line  with  the  data  was 
to  vary  the  values  of  po  for  the  different  phases.  This  is  presumed  to  reflect  the 
effects  of  non-quartz  constituents  of  the  sandstone. 

Above  the  HEL,  the  phase  on  the  sandstone  and  quartz  Hugoniots 
apparently  remains  quartz  until  sufficient  energy  is  available  to  initiate  transition  to 
stishovite.  Although  part  of  the  Hugoniot  is  in  the  coesite  stability  field,  the 
available  data  show  no  evidence  of  coesite  formation  in  the  shock  state.  A  wide 
mixed-phase  (MP)  region  exists  on  the  Hugoniot  between  the  onset  and  completion 
of  stishovite  formation.  The  lower  and  upper  pressure  boundaries  to  this  region 
depend  on  the  value  of  Poo-  The  lower  boundary  on  the  MP  regime  occurs 
considerably  above  the  pressure  at  which  the  Hugoniot  enters  the  stishovite  stability 
field,  but  decreases  as  Poo  decreases,  suggesting  that  the  transition  is  kinetically 
inhibited  (see,  e.g.,  Podurets  and  Trunin  [1987]).  We  can  model  the  fraction  f  of 
the  higher-pressure  phase  by  assuming  an  Arrhenius  rate  constant  for  the  transition 
and  that  the  timescale  for  the  transition  is  relatively  constant  (i.e.,  that  the  thickness 
of  the  shock  front  is  proportional  to  Uj.); 

f  =  Ae-"‘/"^  (12) 

where  the  preexponential  factor  A  and  the  activation  energy  E^  are  determined 
empirically  from  the  experimental  data.  The  assumption  here  is  that  no  transition 
occurs  until  the  phase  boundary  is  crossed,  and  then  is  controlled  by  (12)  after  the 
phase  boundary  is  crossed.  The  issue  of  kinetics  raises  the  possibility  that 
specimen  size  could  affect  the  behavior  of  the  Hugoniot.  However,  this  was 
experimentally  tested  by  Podurets  et  al.  [1976],  who  concluded  that  specimen  size 
does  not  affect  shock  propagation  or  the  position  of  the  P-V  Hugoniot  curve.  For 
the  present  model,  we  assume  that  the  quartz-stishovite  transition  occurs  at  a  set 
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Figure  3.  Equation  of  stale  fits  to  shock  Hugoniot  data  for  quartz  and  fused  quartz  in  the  high  coordination  number  liquid 
stability  region. 


Table  3.  Equation  of  state  properties  of  phases  for  pure  Si02  and  Coconino  Sandstone. 
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pressure  of  12  GPa.  To  simplify  the  calculation  from  a  practical  standpoint,  we 
also  assume  that  the  relevant  temperature  T  is  that  for  the  lower  pressure  phase  at 
the  correct  shock  pressure.  From  the  present  data  and  data  presented  by  Chhabildas 
[1986]  for  quartz  and  Van  Thiel  [1977]  for  Coconino  Sandstone  with  Poo  =  1.98 
Mg/m^,  we  find  that  the  preexponential  factor  has  a  constant  value  of  A  =  7.286, 
while  the  activation  energy  obeys 

E,(kJ/mol)  =  89.36-71.97(po/poo)  (13) 

where  po  is  the  value  for  quartz. 

The  MP  region  is  modeled  by  a  volume  weighted  average  of  the  Hugoniot 
curves  of  the  low  and  high  pressure  phases  in  the  P-V  plane: 

^MP  “  ^HPP  ■*■(!“  ^)^LPP  (1^) 

where  the  subscripts  LPP,  MP,  and  HPP  refer  to  the  low  pressure  phase,  mixed 
phase,  and  high  pressure  phase,  respectively.  This  is  not  thermodynamically 
rigorous,  but  is  convenient  for  the  calculations  and  is  a  very  good  approximation, 
considering  that  thermal  effects  on  V  are  fairly  small  at  high  pressures. 

The  transition  region  between  stishovite  and  the  high-density  liquid  can  be 
treated  in  the  same  way  as  that  between  quartz  and  stishovite.  The  metastable 
Hugoniot  curves  in  this  region  are  very  similar  for  the  cases  where  data  exist  (i.e., 
on  the  quartz  and  fused  quartz  Hugoniot  curves)  so  that  the  transition  region  is  not 
readily  apparent  in  the  P-V  data.  However,  the  transition  region  is  very  apparent  in 
the  P-T  data  of  Lyzenga  et  al.  [1983].  We  use  a  straight  line  fit  to  our  interpreted 
phase  boundary  positions  from  Lyzenga  et  al.  [1983]  as  the  crudest  representation 
of  the  stishovite-liquid  phase  boundary: 

T  =  4348-12.121P  (15) 

for  T  in  K  and  P  in  GPa.  Their  data  show  that,  similar  to  the  quartz-stishovite 
transition,  the  stishovite-liquid  transition  does  not  begin  until  the  Hugoniot  curve  is 
well  into  the  liquid  stability  region.  The  transition  then  occurs  over  a  fairly  narrow 
pressure  range.  The  details  of  the  transition  region  cannot  be  determined  with 
confidence  from  the  existing  data,  so  we  have  chosen  to  express  the  fraction  of 
liquid  by  (12),  with  A  =  6.245  x  10^  and  an  activation  energy  given  by 
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E,(kJ  /  mol)  =  404  +  0.3 12poo 


(16) 


for  Poo  in  kg/m^. 

Figure  4  shows  the  Hugoniot  curves  we  calculate  for  quartz  and  Coconino 
sandstone  using  the  EOS  parameters  in  Table  3  and  the  kinetics  model  discussed 
above. 


3.  Release  Adiabats. 

The  release  adiabats  of  each  of  the  Hugoniot  phases  can  be  calculated  from 
the  EOS  parameters  listed  in  Table  3,  under  the  assumption  that  the  constant  volume 
specific  heat,  Cy,  does  not  vary  significantly,  which  is  a  reasonable  approximation. 
For  a  given  density  p^  on  the  isentrope  anchored  to  a  shock  state  at  density  p^  and 
shock  pressure  the  pressure  P,  on  the  release  isentrope  is 


p,  =  p„  +  (P„  -  Ps)(V„/V,)-‘  expC f  “  (7/V)dV 


(17) 


where  is  the  pressure  in  the  principal  isentrope  corresponding  to  V  =  V,  in  (7)  - 
(9).  The  internal  energy  difference  between  this  and  the  Hugoniot  state  (i.e.,  the 
energy  converted  from  internal  to  kinetic  energy  by  doing  work  to  expand  the 
material  during  the  release  process)  is 

E„-E,  =  E„-E3,-^(P.-P„)  (18) 


where  Eh  is  given  by  (4)  The  internal  energy  not  converted  into  kinetic  energy  upon 
release  is  deposited  as  heat.  This  irreversible  energy  deposition  causes  the  shock 
wave-release  wave  couple  to  decay  as  energy  required  to  drive  the  wave 
propagation  is  sapped  away. 

The  particle  velocity  of  material  released  from  the  shocked  state  is  related  to 
the  P-V  release  path  via  the  Riemann  integral  formula  [Rice  et  al.,  1958]: 
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Hgure  4.  Model  Hugoniots  and  release  adiabats  for  (a)  quartz  and  (b)  Coconino 
Sandstone. 
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where  the  subscript  rel  refers  to  the  isen tropic  release  path.  An  additional 
contribution  to  u^  can  come  from  the  energies  of  transition  associated  with  phase 
changes  during  release.  We  have  assumed  that  the  total  internal  energy  is  constant 
with  latent  heat  affecting  sensible  heat  rather  than  macroscopic  kinetic  energy.  The 
simplest  model  assumes  that  the  relative  fractions  of  the  low  and  high  pressure 
phases  are  kinetically  frozen  at  their  Hugoniot  state  values.  The  available  evidence 
from  recovered  shocked  quartz  and  from  the  present  data  is  that  the  stishovite 
transforms  to  a  diaplectic  glass  with  an  STP  density  of  2.27  Mg/m^  and  that  dense 
amorphous  Si02  transforms  to  the  more  familiar  low-pressure  glass  (p  =  2.204 

Mg/m^)  or  liquid  at  low  pressures.  Based  on  the  available  data,  we  assume  this 
happens  at  a  critical  pressure  of  ~5.1  GPa  for  both  phases.  In  both  cases,  we 
have  assumed  that  the  path  below  P,  is  a  straight  line  in  P-V  space  from  the  volume 
at  Pj  to  the  STP  volume  of  fused  quartz.  This  ignores  the  effects  of  thermal 
expansion,  but  the  error  introduced  is  insignificant.  As  with  the  Hugoniot,  we 
model  the  release  points  on  the  mixed-phase  release  paths  by  assuming  a  mass- 
weighted  average  of  the  volumes  of  the  low  and  high  pressure  phases  via  (14). 
Figure  4  presents  release  adiabats  calculated  with  this  model. 


C.  Limestone. 

1.  Experimental  Data. 

Table  4  presents  the  data  obtained  for  Solenhofen  and  Bedford  limestones  in 
this  study.  Figures  5  and  6  present  these  data  along  with  data  from  Ahrens  and 
Gregson  [1964],  Van  Thiel  [1977],  and  Tyburczy  and  Ahrens  [1986]. 

The  low  pressure  behavior  of  CaC03  is  very  complicated.  At  least  four 

phases  are  known  which  seem  to  be  only  slightly  different  in  internal  energy.  In 
addition  to  calcite  I  and  aragonite,  at  least  two  phases,  calcite  D  and  calcite  HI,  are 
observed  in  static  compression,  with  the  calcite  I-II  and  II-III  transitions  occurring 
at  -1.5  GPa  and  2.2  GPa,  respectively,  at  298  K  [Fiquet  et  al.,  in  press].  Static 
compression  data  of  Adadurov  et  al.  [1960]  also  show  a  very  slight  volume 
discontinuity  at  -1.16  GPa,  although  this  could  be  due  to  some  experimental  effect 
Complicating  this  issue  is  that,  with  the  exception  of  calcite  I,  different  investigators 
disagree  markedly  in  their  determinations  of  the  properties  of  these  phases. 

Low  pressure  shock  wave  data  for  calcite  and  calcite  rocks,  including  the 
data  obtained  in  the  present  study,  show  not  only  an  elastic  precursor  wave  at  low 
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Table  4.  Shock  wave  and  release  adiabat  data  for  Solenhofen  and  Bedford  limestones. 
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Table  4.  (continued) 
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Figure  5.  Shock  Hugoniot  and  release  adiabat  data  and  Hugoniot  based  on 
empirical  U»-Up  fits  for  Solenhofen  Limestone,  (a)  U,-Up  plane,  (b)  P-V  plane. 
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Figure  6.  Shock  Hugoniot  and  release  adiabat  data  and  Hugoniot  based  on 
empirical  U,-Up  fits  for  limestone  from  Bedford,  Illinois,  (a)  Uj-Up  plane,  (b)  P-V 
plane. 
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shock  stresses,  but  successive  formation  of  at  least  two  and  perhaps  three  high- 
pressure  phases,  detected  by  the  sequential  arrival  of  multiple  shock  fronts. 

Although  the  data  tend  to  be  somewhat  scattered,  the  stress  levels  of  these  different 
waves  in  the  multiple-wave  experiments  are  similar  to  those  characterizing  calcite  n 
and  ni  formation  in  static  compression  experiments. 

As  with  sandstone,  the  velocity  of  the  elastic  precursor  wave  is  not  strongly 
correlated  with  particle  velocity  and  is,  to  first  order,  a  function  only  of  the  initial 
density  of  the  limestone.  We  have  combined  the  data  of  Ahrens  and  Gregson 
[1964]  for  the  first  wave  velocities  in  polycrystalline  calcite  rocks  with  those  from 
our  experiments  to  obtain  a  data  set  for  rocks  with  poo  in  the  range  from  2.327 
Mg/m^  to  2.619  Mg/m\  A  fit  of  Uji  as  a  function  of  poo  gives 

U„  =-80.256 -h63.217poo  -11.617poo'  (20) 

where  U^i  is  in  km/s  and  poo  is  in  Mg/m^. 

2.  Equation  of  State. 

Because  of  the  confused  nature  of  the  CaC03  system  at  low  pressures,  we 
can  provide  only  a  very  approximate  treatment  in  that  region.  One  very  important 
complication  is  that  the  final  shock  states  of  calcite  rocks  below  ~25  GPa  cannot  be 
matched  by  any  published  equation  of  state  for  any  single  phase  of  CaC03.  We 

have  chosen  to  model  the  lowest  pressure  portions  of  the  Hugoniots  of  low- 
porosity  calcite  rocks  as  a  mixture  of  calcite  I  and  calcite  III,  where  we  identify 
calcite  III  as  the  phase  observed  by  both  Piquet  et  al.  [1993]  and  Vizgirda  and 
Ahrens  [1982].  In  the  case  of  Bedford  limestone,  (porosity  ~  12%),  aragonite  is 
used  instead  of  calcite  III,  as  explained  below. 

The  identity  of  the  phase  comprising  the  final  shock  state  just  above  this 
mixed  phase  region  is  uncertain.  The  shock  compression  data  for  single  crystal 
calcite  from  Ahrens  and  Gregson  [1964]  seem  to  indicate  that  the  final  phase  is  not 
calcite  ID,  but  a  higher  pressure  phase  called  calcite  IV  by  Tyburczy  and  Ahrens 
[1986].  However,  Piquet  et  al.  [1993]  present  new  X-ray  diffraction  data  which 
indicate  that  calcite  III  has  the  same  properties  as  calcite  IV  given  by  Vizgirda  and 
Ahrens  [1982].  The  previous  EOS  parameters  for  calcite  III  are  called  into  question 
because  the  phase  coexistence  regions  noted  by  Piquet  et  al.  [1993]  can  introduce 
errors  into  non-X-ray  volume  determinations  used  in  the  earlier  studies.  Piquet  et 
al.  [1993]  do  note  the  appearance  of  new  X-ray  lines  and  a  possible  calcite  IV  phase 
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at  about  the  same  pressure  as  found  by  Ahrens  and  Gregson  [1964],  but  the 
properties  of  this  phase  under  shock  compression  seem  to  be  consistent  with  calcite 
TTT,  so  the  situation  is  at  present  confused.  A  further  complication  is  that  the 
Bedford  Limestone  Hugoniot  above  this  point  is  more  consistent  with  the  properties 
of  aragonite  than  calcite  in.  The  issue  will  likely  require  further  careful  static  X-ray 
diffraction  experiments  at  elevated  temperature  and  pressure  to  map  out  the  phase 
diagram  before  it  can  be  resolved. 

Based  on  the  available  data  and  the  EOS  parameters  presented  by  Piquet  et 
al.  [1993]  and  Vizgirda  and  Ahrens  [1982],  we  assume  shock  states  in  the  25-40 
GPa  range  for  calcite  to  be  a  single  phase  which  we  identify  as  calcite  HI.  Based  on 
the  published  data  and  previous  determinations  of  EOS  parameters,  we  propose  that 
calcite  III  is  described  by  Poo  =3.100  Mg/m^  Kgo  =  96.5  GPa,  K’  =  4.0,  Yo  = 
1.183,  n  =  0.6,  and  E^,  =  0.02  MJ/kg.  These  parameters  also  describe  the 
equivalent  portions  of  the  aragonite  and  relatively  low  porosity  limestone  Hugoniots 
reasonably  well.  For  limestones  with  porosities  as  much  as  15%,  the 
corresponding  region  is  better  fit  by  aragonite. 

Above  the  calcite  HI  and  aragonite  single  phase  regions,  the  Hugoniot  state 
appears  to  change  to  yet  another  phase.  We  fit  the  shock  compression  data  for  P  > 
100  GPa  for  calcite  and  limestones  to  obtain  =  3.377  Mg/m 3,  Kjo  =137  GPa, 

K'  =  3.4,  Yo  =  2.5,  n  =  1.0,  and  E^  =  2.142  MJ/kg.  We  presume  that  this  phase 
probably  is  a  dense  melt,  given  the  large  STP  transition  energy  from  calcite  I. 

Table  5  summarizes  the  EOS  parameters  we  have  adopted  for  each  phase 
considered  in  this  study.  We  assumed  in  all  cases  that  Cv  =  3k/atom.  Calcite  n  is 
not  actually  used  in  our  study,  but  is  presented  for  completeness.  Another 
complication  with  calcite  rocks  is  that  no  combination  of  the  phases  described  in 
Table  5  will  give  the  correct  Hugoniots  for  highly  porous  rocks  such  as  chalk  (see, 
e.g.,  Kalashnikov  et  al.  [1973],  and  Tyburczy  and  Ahrens  [1986]).  Because  of  the 
large  internal  energies  induced  by  shock  waves  in  porous  materials,  it  is  probable 
that  the  phases  encountered  on  the  Hugoniots  of  these  materials  are  not  the  same  as 
those  found  on  the  Hugoniots  of  the  denser  rocks  studied  here.  We  attempted  to  fit 
separate  equations  of  state  for  the  data  for  chalk,  but  still  could  not  get  consistent 
results.  We  suggest  that  the  Hugoniot  states  of  very  porous  calcite  rocks  may  be 
affected  by  decomposition  of  the  CaC03.  Such  behavior  will  exacerbate  the 
difficulties  of  developing  accurate  models  of  the  behavior  of  porous  CaC03  rocks. 
At  present,  we  advise  that  the  models  developed  in  this  study  be  restricted  to 
limestones  with  <20%  porosity. 
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Table  5.  Equation  of  state  properties  of  CaC03  phases. 


Property 

Calcite  I 

Calcite  11 

Aragonite 

Calcite  HI 

Dense  Liquid 

po  (Mg/m^) 

2712 

2522 

2931 

3100 

3377 

Kso  (GPa) 

74.405 

13 

66.73 

96.5 

139 

K’ 

5.42 

4 

4 

4 

3.4 

yo 

0.554 

1.5 

1.80 

1.183 

2.5 

n 

1 

1 

1 

0.5 

1 

E^(MJ/kg) 

0 

0 

-0.0006 

0.02 

2.142 

Cv  (J/kg/K) 

1246 

1246 

1246 

1246 

1246 

Source 

(1) 

(2) 

(3) 

(4) 

(5) 

(1)  Sumino  and  Anderson  [1984] 

(2)  Fit  to  data  from  Piquet  et  al.  [1993] 

(3)  Calculated  from  Salje  and  Viswanathan  [1976]  and  Robie  et  al.  [1978] 

(4)  Piquet  et  al.  [1993]  and  Vizgirda  and  Ahrens  [1982] 

(5)  This  study 
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As  with  the  Si02  system,  we  model  the  mixed-phase  regions  of  the 
Hugoniot  by  assuming  kinetically-inhibited  transformations  on  the  equilibrium 
phase  boundary  is  crossed.  The  primary  reason  for  taking  this  assumption  is  the 
fact  that  formation  of  the  high  pressure  phase  identified  as  calcite  in  by  Piquet  et  al. 
[1993]  above  2.2  GPa  is  not  complete  until  ~25  GPa  for  single-crystal  calcite  under 
shock  compression.  Using  the  EOS  parameters  in  Table  5  to  calculate  the  Hugoniot 
temperatures  corresponding  to  the  P-V  data,  we  fit  Arrhenius  parameters  for  (12)  to 
get  A  =  5.848  and 

E3  =  62.713-20.651poo  (21) 

for  Pqq  in  Mg/m^  and  E^  in  kJ/mol  for  the  conversion  of  calcite  I  to  calcite  HI  or 
aragonite.  For  the  calcite  III  -  liquid  and  aragonite  -  liquid  mixed  phase  regions,  we 
get  A  =  2.919  and 

E3  =  382. 453  - 1 30. 448poo  (22) 

for  Poo  in  Mg/m^  and  E^  in  kJ/mol.  For  the  calcite  I-aragonite  and  calcite  I-calcite  HI 
phase  boundaries,  we  take  the  overall  trends  of  the  data  presented  by  Jamieson 
[1957]: 

Ti_^=164.7P-3.6  (23) 

Tj.iji  =78.0P-1.3  (24) 

for  T  in  K  and  P  in  GPa.  For  the  provisional  calcite  Ill-liquid  and  aragonite-liquid 
phase  boundaries,  we  get 

TAr-uq=745-H48P  (25) 

Tin.u<,=401-h25.8P  (26) 

for  T  in  K  and  P  in  GPa.  These  phase  boundaries  are  extremely  uncertain,  being 
based  solely  on  the  behavior  of  the  P-V  Hugoniots  and  the  temperatures  calculated 
from  the  EOS  parameters  given  in  Table  5.  Figure  7  presents  the  Hugoniots  of 
Solenhofen  and  Bedford  limestones  calculated  with  this  model. 
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Volume  (mV^^9) 


Figure  7.  Model  Hugoniots  and  release  adiabats  for  (a)  Solenhofen  and  (b) 
Bedford  limestones. 
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3.  Release  behavior. 

Measured  points  on  release  adiabats  from  the  calcite  III  portion  of  the 
Hugoniot  seem  to  require  reversion  to  calcite  I  with  energy  approximately  equal  to 
the  STP  transition  energy  being  converted  to  kinetic  energy.  One  problem  with  this 
interpretation,  however,  is  that  the  phase  reversion  takes  place  at  a  pressure  higher 
than  would  be  expected  based  solely  on  the  CaC03  phase  transition  pressures.  We 
may  be  detecting  partial  decomposition  of  calcite  III  or  formation  of  a  phase  not  yet 
identified.  TTie  final  volumes  of  samples  released  from  the  calcite  HI  region  to  P  = 

0  are  large  and  suggest  that  CO2  gas  is  being  evolved  at  very  low  pressures.  There 

is  no  evidence  for  this  behavior  at  P  =  1.41  GPa,  so  that  the  decomposition 
probably  occurs  at  P  <  1  GPa. 

In  contrast  with  release  from  calcite  III,  the  liquid  phase  seems  to  be 
kinetically  inhibited  from  freezing  down  to  release  state  pressures  as  low  as  4  GPa 
and  follows  the  release  isentropes  described  by  (16).  We  suspect  that  this  phase 
also  undergoes  decomposition  with  evolution  of  CO2  at  P  <  1  GPa. 

As  with  Si02,  we  model  release  from  mixed-phase  regions  of  the  Hugoniot 
as  release  from  the  meiastable  Hugoniots  of  the  two  end-member  phases  at  the  same 
pressure.  For  this  puipose,  we  assume  that  release  from  the  calcite  I  metastable 
Hugoniot  is  the  simple  release  isentrope  described  by  (16).  Release  adiabats 
described  by  this  model  are  shown  in  Figure  7. 

D.  Calcareous  Shale 
1.  Samples. 

A  calcareous  shale  was  collected,  by  us,  for  this  study.  It  was  obtained 
from  the  lower  Chambless  Limestone  (Cambrian)  in  the  southernmost  Marble 
Mountains  of  the  central  Mojave  Desert  in  California.  Volatile-free  oxide 
composition  was  determined  via  inductively-coupled  plasma  mass  spectroscopy  and 
X-ray  fluorescence  spectrometry.  CO2  in  carbonates  was  determined  by  reaction 

with  acid  and  volumetric  measurement  of  gas  evolved.  The  composition  is  given  in 
Tables  6  and  7.  The  samples  had  a  mean  nonporous  (crystal)  density  of  2.615 
Mg/m^and  a  mean  bulk  density  of  2.576  Mg/ml  Ultrasonic  measurements  of 
longitudinal  and  transverse  wave  velocities  give  Vp  =  4.505  km/s  and  V,  =  2.753 
km/s,  respectively. 


26 


Table  6. 

Oxide  composition  analyses  of  Chambless  calcareous  shale  from  Marble  Mountains, 
central  Mojave  Desert,  California. 


Oxide 

Weight  %  (ICP)a 

Weight  %  (XRF)b 

SiOa 

44.59 

44.20 

TiOj 

0.84 

0.76 

A1203 

9.68 

9.42 

FejOs 

4.77 

4.73 

CaO 

19.09 

19.30 

MgO 

1.33 

1.39 

Na20 

0.33 

0.13 

K2O 

2.04 

2.25 

C02= 

14.34 

14.34 

H20d 

4.21 

3.46 

Total 

100.94 

99.98 

(a)  Indued vely-couple  plasma  source  emission  spectrometry 

(b)  X-ray  fluorescence  spectrometry 

(c)  Determined  by  gas  evoludon  on  treatment  with  acid. 

(d)  Determined  by  thermal  desorption  and  decomposidon  (different  samples). 
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Table  7 

Mineral  composition  of  calcareous  shale  used  in  this  study,  estimated  from  oxide 
composition  and  X-ray  diffraction  analysis. 

Mineral  Weight  Percent 


Quartz  32.1 

Calcite  32.6 

Clays  and  Micas  34.4 

Zeolitic  Water  1 . 1 

Total  100.2 
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2.  Calcareous  Shale,  Experimental  Results 

Table  8  and  Figure  8  present  our  results  for  calcareous  shale.  The  U^-Up 
data  can  be  broken  down  into  three  regions  and  fit  by  equation  (1)  truncated  to  first 
order.  The  fit  parameters  are  Q  =  3.803±0.097  km/s  and  s  =  1.463±0.1 15  for  Up 
<  1.12  km/s,  Co  =  4.853±0.1 15  km/s  and  s  =  0.523±0.073  for  1.12  <Up<  1.98 
km/s,  and  Cq  =  2.405±0.139  km/s  and  s  =  1.761±0.054  for  Up  >  1.98  km/s.  The 
high  compressibility  shown  at  low  pressures  in  the  P-V  data  and  the  definitely 
defined  segments  of  the  U^-Up  Hugoniot  suggest  that  the  lowest  pressure  data 
represent  Hugoniot  states  entirely  within  the  low-pressure  phase  of  the  shale  and 
the  highest  pressure  data  represent  states  entirely  within  a  high  pressure  phase.  The 
intervening  region  represents  a  mixed-phase  regime  where  the  formation  of  the  high 
pressure  phase  is  incomplete,  so  that  the  Hugoniot  states  fall  on  a  mixing  line 
between  the  metaslable  extensions  of  the  low  and  high  pressure  phase  Hugoniot 
curves.  We  should  note  here  that  the  term  "phase"  is  not  strictly  correct,  since  the 
rock  is  a  mixture  of  several  different  mineral  phases.  Rather,  the  term  "low 
pressure  phase"  in  the  present  context  is  used  to  refer  to  a  material  composed  of  the 
same  mineral  phases  as  the  original  sample,  whereas  the  term  "high  pressure  phase" 
refers  to  a  material  in  which  some  or  all  of  the  constituent  minerals  have 
transformed  to  high  pressure  polymorphs.  The  term  "mixed  phase"  in  the  present 
paper  refers  to  a  material  in  which  minerals  have  partially  transformed  to  high 
pressure  polymorphs,  but  the  transformation  is  incomplete. 

In  the  lowest  pressure  experiments,  the  plastic  shock  wave  in  the  shale  was 
preceded  by  an  elastic  precursor.  If  we  assume  that  the  velocity  i/«(  of  elastic 
precursor  should  be  identical  to  Vp  under  infinitesimal  strain,  then  the  elastic 
precursor  obeys  t/,,  =  4. 505  -I-  L  2Up  for  and  Up  in  km/s. 
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Particle  Velocity  (km/s) 


Volume  (mV^Q) 


Figure  8.  Experimental  data  for  calcareous  shale:  (a)  U^-Up  projection  with 
empirical  fits  to  the  data;  (b)  P-V  projection  of  data  and  Hugoniot  calculated  from 
U,-Up  fit. 
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Tables 

Experimental  E>ata  for  Calcareous  Shale. 
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3.  Theoretical  Equation  of  State  of  Shale. 

The  shale  we  studied  is  relatively  carbonate  rich  (32.6%  CaCOj  by  weight). 
The  carbonate  content  of  shales  can  be  quite  variable,  so  we  have  concentrated  here 
on  a  carbonate-rich  shale  and  use  published  data  for  carbonate-poor  shales  to  obtain 
two  different  equations  of  state.  We  used  the  composition  (Table  6),  on  an  HjO 
free  basis,  to  calculate  the  equivalent  mineral  mixture  model  of  the  rock,  and  then 
accounted  for  H2O  in  the  analysis  and  corundum  in  the  norm  by  converting 
feldspars  and  other  minerals  to  clays  and  micas.  We  represent  H2O  bound  in  the 
structure  of  clays  and  micas  as  a  separate  constituent.  Table  7  gives  the  resulting 
approximate  mineral  composition.  The  rock  is  composed  of  roughly  equal  amounts 
of  quartz,  clays,  and  calcite.  This  result  is  confinned  by  powder  X-ray  diffraction 
analysis.  The  mean  Archimedian  value  for  po  of  the  low  pressure  phase  is  2.615 
Mg/m^  with  calculated  values  of  yb  =  0.73,  n  =  1.52,  Cv=  1396  Jkg  ’K-'.  We 
held  Kos  at  30.9  GPa,  calculated  from  (12),  and  find  that  the  fit  to  the  low  pressure 
data  not  affected  by  an  elastic  precursor  gives  fC  =  8.15.  The  calculated  parameters 
for  the  high  pressure  phase  are  /o  =  1-2,  n  =  1.49,  £,r=  .42  MJ/kg,  and  Cv=  1396 
Jkg-'K-‘.  A  fit  to  the  four  highest  pressure  shock  compression  data  gives  Po  = 
3.659  Mg/m^  Kos=  265  GPa  and  K'  =  3. 14. 

To  study  relatively  carbonate-poor  shales,  we  have  combined  the  data  of 
Trunin  et  al.  [1988]  and  Al'tshuler  and  Pavlovskii  [1971]  for  Russian  shales  and  of 
Olinger  [1978]  for  kerogen-free  Chemung  Shale.  Trunin  et  al.  [1988]  give  a 
chemical  analysis  of  their  shale,  which  we  use  as  representative  of  all  the  samples, 
since  the  three  data  sets  are  the  same  within  the  scatter  of  the  data.  The  analysis 
give  ranges  of  abundances  of  different  constituents,  so  that  our  estimation  of  the 
mineral  content  of  the  shales  is  only  approximate.  They  also  do  not  report  an  H2O 
analysis  but,  based  on  the  normative  comndum  abundance  and  the  densities  of  the 
rocks,  the  dominant  phyllosilicates  are  probably  micas,  rather  than  clays  with  the 
rock  containing  only  -8%  CaCOs.  Table  9  gives  the  estimated  mineral  composition 
of  these  shales.  We  estimated  po  for  the  low  pressure  phase  on  the  assumption  that 
all  the  normative  corundum  was  in  muscovite  and  get  po  =  2.813  Mg/m^.  Other 
calculated  parameters  are  %  =  0.79,  n  =  1.73,  and  Cv=  1 159  Jkg  'K  *.  To  avoid 
possible  problems  caused  by  elastic  precursors  in  the  original  data,  we  fit  the  low 
pressure  phase  EOS  to  the  data  in  the  pressure  range  from  10  GPa  to  16  GPa  and 
get  Xos=  61.1  GPa  and  K  =  4.52.  Calculated  parameters  for  the  high  pressure 
phase  are  /o  =  0.93,  n  =  1.92,  Cv=l  159  Jkg"‘K‘’ ,  and  E,r  =  0.62  MJ/kg.  A  fit  to 
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data  in  the  pressure  range  from  55  GPa  to  120  GPa  gives  po  =  3.740  Mg/m^  Kos= 
184  GPa  and  /T  =  3.24. 
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Table  9 

Estimated  mineral  composition  of  carbonate-poor  shale  studied  by  Trunin  et  al.  [1988]. 


Mineral 

Weight  Percent 

Quartz 

40 

Plagioclase 

15 

Pyroxene 

13 

Calcite 

8 

Muscovite 

20 

Magnetite 

4 

Total 

100 
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4.  Release  Adiabats  of  Shale 

Figure  9  presents  the  experimental  P-V  release  adiabat  points  obtained  in 
this  study.  In  all  cases,  the  release  adiabats  from  the  high  pressure  region  of  the 
Hugoniot  parallel  the  P-VHugoniot  at  high  pressures,  with  deviation  to  larger 
volumes  at  lower  pressures.  Release  from  the  mixed-phase  regions  seems 
generally  to  result  in  large  volume  increases  probably  related  to  reversion  to  the  low 
pressure  phase.  Calcareous  shale  shows  a  large  volume  increase  upon  release, 
especially  in  the  mixed-phase  region,  similar  to  the  carbonate  rocks  and  minerals 
[Vizgirda  and  Ahrens,  1982;  Anderson  and  Ahrens,  1995].  The  large  volume 
increase  at  low  pressures  strongly  suggest  that  the  high  pressure  phase  of  shale 
reverts  to  the  low  pressure  phase  or  to  some  other  large-volume  phase  upon  release. 

As  with  the  sandstone  and  limestone,  the  experimental  release  points  can  be 
used  to  check  the  validity  of  our  equations  of  state.  The  points  shown  are  actually 
lower  limits  on  V,  since  the  straight-line  approximation  of  Lyzenga  and  Ahrens 
[1978]  produces  the  lower  extremum  of  thermodynamically  allowable  volume 
increase  upon  adiabatic  release  with  no  phase  change.  For  a  rigorous  comparison, 
one  would  need  to  calculate  the  particle  velocity  increase  for  a  test  release  path  using 
the  Riemann  integral  and  compare  that  to  the  observed  particle  velocity  increases. 
However,  since  the  release  state  particle  velocities  often  have  fairly  large 
uncertainties,  the  straight-line  approximation  of  the  release  path  is  sufficient  for  the 
present  study.  We  again  make  the  assumption  that  Cy  is  constant  so  that  eq.  (17) 
can  be  used  to  calculate  theoretical  release  adiabats.  Figure  10  shows  the  release 
adiabats,  calculated  under  the  assumption  that  the  high  pressure  phase  is  "frozen  in" 
at  the  shock  state  mass  fraction  by  kinetics.  Reversion  of  the  high  pressure  phase 
to  some  other  lower  density  phase  is  manifested  in  the  figures  by  experimental 
release  paths  at  lower  pressures  that  are  shallower  than  the  calculated  paths. 
However,  the  higher  pressure  data  generally  show  good  agreement  with  the 
calculated  release  adiabats. 
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Figure  9.  Calculated  Hugoniot  and  release  adiabats  for  calcareous  shale,  compared 
with  the  experimental  data. 


37 


Volume 


OSL  001  09  0 

(Ddo)  9jnssajc| 

Figure  10.  Calculated  Hugoniot  and  release  adiabats  for  carbonate-poor  shale, 

compared  with  the  experimental  data. 
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III.  DISCUSSION 


Now  that  a  model  for  the  total  shock-release  cycle  has  been  developed  we 
can  apply  it  to  calculate  the  fraction  f  of  the  shock  internal  energy  Ejj  that  is 
irreversibly  deposited  in  the  material  for  a  given  shock  pressure  using  the  equation 


f  = 


Eh-E. 

E„ 


(27) 


where  Eh  and  the  quantity  (Eh-Ej)  are  given  by  (4)  and  (18),  respectively.  We  have 
assumed  here  that  the  effects  of  the  elastic  precursor  at  low  pressures  can  be 
ignored.  The  resulting  estimates  for  f  are  presented  in  Figures  1 1  and  12.  For 
Coconino  Sandstone,  f  approaches  1  at  low  shock  pressures  because  of  pore 
collapse,  but  rapidly  drops  to  -0.44  at  18  GPa,  rising  rapidly  to  -0.65  at  60  GPa 
and  varying  slowly  thereafter,  decreasing  to  -0.59  at  200  GPa.  In  the  case  of  the 
limestones,  we  have  ignored  the  rapid  expansion  due  to  CO2  gas  evolution  at 
extremely  low  pressures,  so  the  values  of  f  given  are  upper  limits  in  the  case  of 
release  to  P=0.  At  low  pressures,  as  with  sandstone,  f  approaches  1  due  to  pore 
collapse,  but  drops  rapidly,  the  rate  of  the  drop  being  dependent  on  the  initial 
density.  Because  of  differences  in  the  phases  on  the  Hugoniot  (aragonite  in 
Bedford  limestone  versus  calcite  HI  in  Solenhofen  Limestone)  and  the  difference  in 
Poo,  the  P-f  curves  behave  very  differently  in  detail  below  120  GPa.  Above  120 
GPa,  both  curves  are  very  similar,  but  with  slightly  greater  energy  deposition 
(higher  f )  in  the  more  porous  Bedford  limestone,  as  should  be  expected. 

There  are  some  fundamental  observations  that  should  be  made  about  the 
curves  in  Figure  11.  The  detailed  structure  of  the  curves  should  be  taken  with  a  fair 
amount  of  skepticism,  since  the  models  for  the  phase  transition  kinetics  and  the 
phase  boundaries  are  based  on  sparse  data.  However,  the  relative  positions  of  the 
curves  is  a  property  only  of  the  equations  of  state  thus  should  be  a  robust  result. 
Since  a  higher  value  of  f  is  positively  correlated  with  the  irreversible  energy 
deposition  and  thus  the  rate  of  decay  of  the  shock  wave,  rocks  characterized  by 
high  f  values  will  tend  to  decouple  explosions  from  far  field  observables  more 
efficiently  than  would  rocks  with  low  f  values.  In  particular,  the  conclusion  from 
the  curves  presented  in  Figure  1 1  is  that  porous  silicate  rocks  such  as  sandstones 
are  much  better  suited  to  decoupling  that  are  carbonate  rocks.  This  is  aside  from 
other 
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Figure  1 1 .  Fraction  f  of  Rankine-Hugoniot  internal  energy  which  is  irreversibly 
deposited  as  heat  in  Coconino  Sandstone  and  Solenhofen  and  Bedford  limestones, 
shown  as  a  function  of  peak  shock  pressure. 
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Figure  12.  Fraction  f  of  Hugoniot  internal  energy  irreversibly  deposited  in  shale 
and  slate,  calculated  assuming  retention  of  the  high  pressure  phase  upon  release. 
Actual  energy  deposition  will  be  smaller. 
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issues,  such  as  the  need  to  contain  gas  evolved  from  decomposition  of  CaCOs  as 
release  pressures  approach  zero.  Figure  12  summarizes  the  values  of  f  for  the 
shales  and  slates. 

An  important  effect  which  has  been  ignored  here  is  the  effect  of  the  elastic 
precursor  wave  that  precedes  the  plastic  shock  wave  at  low  pressures.  The  actual 
values  of  Eh  for  our  limestone  experiments  which  showed  elastic  precursors  and 
intermediate  waves  were  ~35%  and  ~65%  for  Bedford  and  Solenhofen  limestones, 
respectively,  of  the  value  that  is  calculated  for  the  same  final  particle  velocity  in  a 
purely  plastic  shock  such  as  we  have  used  in  the  present  calculations.  For  the  data 
from  Ahrens  and  Gregson  [1964]  for  Coconino  Sandstone,  the  actual  internal 
energy  is  -55%  of  the  calculated  plastic  shock  energy.  Also,  material  strength  in 
mitigating  pore  collapse  at  low  pressures,  so  that  the  energy  deposition  is  much 
lower.  Since  the  release  behavior  should  be  the  same,  however,  the  result  is  that  at 
very  low  pressures,  even  in  porous  materials,  f -0.  There  is  ample  experimental 
evidence  that  significant  energy  deposition  is  occurring  at  low  shock  pressures, 
however,  so  that  processes  such  as  fracture  fomiation  must  be  important  to  the 
overall  question  of  coupling  of  the  near-source  behavior  to  far-field  observables. 
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IV.  CONCLUSIONS 


We  have  performed  a  series  of  shock  Hugoniot  and  release  adiabat 
measurements  on  sandstone,  limestone,  and  shale  and  combined  the  results  with 
previous  data  to  develop  models  for  the  shock  and  release  behavior  of  these  rocks. 
The  single-phase  portions  of  the  sandstone  Hugoniot  curve  can  be  calculated  from 
the  equations  of  state  of  quartz,  stishovite,  and  a  high  coordination  number  liquid. 
The  properties  of  quartz  and  stishovite  have  been  extensively  studied  in  the  past. 
We  present  a  new  fit  of  available  data  for  the  dense  liquid  phase  to  get  poo  =  4.150 
Mg/m5,  Kso  =  313.5  GPa,  K'  =  1.8484,  y  =  1.4(V/Vo)'^,  and  £„  =  2.4  MJ/kg.  In 
the  case  of  Coconino  Sandstone,  the  presence  of  a  few  percent  of  non-quartz 
constituents  requires  that  the  STP  densities  of  tlie  relevant  phases  be  reduced  by 
2.4%  to  fit  the  data.  The  intervening  mixed-phases  region  can  be  modeled  as  a 
simple  mixture  of  these  two  relevant  end  member  phases  under  the  assumption  that 
formation  of  the  higher  pressure  phase  is  kinetically  inhibited.  The  fraction  f  of  the 
higher  pressure  phase  is  described  by  f  =  Ae"^'^'^'^.  For  formation  of  stishovite 
from  quartz,  A  =  7.286  and  E3(kJ  /  mol)  =  89.36 -  71. 97(po /poo)  where  poo  and 

Po  are  the  actual  and  theoretical  nonporous  densities  of  the  initial  rock  sample.  For 
formation  of  the  dense  liquid  from  stishovite,  A  =  6.245  x  10^  and  Ea  (kJ/mol)  = 
404  +  0.312poo.  No  formation  of  the  higher  pressure  phase  occurs  until  the  phase 
boundary  is  reached.  As  crude  approximations  we  describe  the  quartz-stishovite 
phase  boundary  by  P  =  12  GPa  and  the  stishoviie-dense  liquid  phase  boundary  by 
T  =  4348-12.121P  for  T  in  K  and  P  in  GPa. 

Release  adiabats  from  the  single  phase  regions  of  the  sandstone  Hugoniot 
are  modeled  by  the  isentropes  of  the  relevant  phases,  as  described  by  the  equation 
of  state.  At  ~5  GPa,  stishovite  begins  to  fonn  a  diaplectic  glass.  We  model  the 
release  adiabat  below  this  pressure  as  a  straight  line  in  P-V  space,  from  the 
stishovite  release  isentrope  at  5  GPa  to  the  STP  diaplectic  glass  density  of  2.27 
Mg/m\  The  dense  liquid  undergoes  a  similar  change  beginning  at  -5  GPa  with 
formation  of  normal  silica  glass  with  an  STP  density  of  2.204  Mg/m^.  The  quartz 
release  adiabat  follows  the  quartz  isentrope  to  P=0.  As  with  the  Hugoniot,  we 
model  the  release  from  the  mixed-phase  regions  of  the  Hugoniot  as  a  mixture  of  the 
relevant  phases,  with  the  proportions  frozen  at  the  Hugoniot  values. 

The  limestone  Hugoniot  curve  is  complicated  by  several  metastable  phases 
and  multiple-wave  shock  structure  at  low  pressures.  Single  phase  portions  of  the 
Hugoniots  of  Solenhofen  and  Bedford  limestones  are  modeled  as  calcite  I  at  low 
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pressures  and  a  dense  liquid  at  high  pressures.  We  have  developed  an  equation  of 
state  for  the  dense  liquid  phase,  with  po  =  3.377  Mg/m^  Kso  =139  GPa,  K’  =  3.4, 
y  =  2.5(VA^o)>  and  =  2.142  MJ/kg.  At  intermediate  pressures,  dense  limestones 
such  as  Solenhofen  seem  to  form  calcite  HI,  while  the  more  porous  Bedford 
limestone  apparently  forms  aragonite.  The  mixed  phase  regions  of  the  Hugoniot 
are  modeled  in  a  manner  analogous  to  that  u.sed  for  the  sandstone  Hugoniot.  The 
parameters  describing  the  fraction  f  of  the  higher  pressure  phase  in  the  mixed  phase 
regions  are  A  =  5.848  and  Ea  =  62.7 13-20.65  Ipoo  for  formation  of  calcite  in  or 
aragonite  from  calcite  I  and  A  =  2.919  and  Ea  =  382.453- 130.448poo  for  formation 
of  the  dense  liquid  from  either  aragonite  or  calcite  111,  for  Ea  in  kJ/mol  and  poo  in 
Mg/m^.  The  phase  boundaries  in  our  model  are  approximated  by  Tj.a,  =  164.7P- 
3.6,  Ti.ni  =  78.00-1.3,  T^^-uq  =  745+48P,  and  Tm.uq  =  401+25.8P. 

As  with  sandstone,  we  model  the  release  adiabats  with  the  isentropes  of  the 
relevant  phases,  mixing  them  at  the  appropriate  ratios  for  release  from  the  mixed- 
phase  portions  of  the  Hugoniot.  There  is  no  evidence  of  formation  of  a  glass  at 
very  low  pressures,  but  release  of  phases  other  than  calcite  I  is  apparently 
accompanied  by  decomposition  with  evolution  of  CO2  gas  at  pressures  below  1 
GPa.  We  have  ignored  this  very  low  pressure  effect  in  our  calculations. 

Equation  of  state  parameters  for  shale  are  calculated  from  the  mineral 
composition,  with  the  bulk  modulus  and  its  pressure  derivative  constrained  by  the 
shock  wave  data.  The  shale,  which  is  composed  of  roughly  equal  amounts  of 
quartz,  calcite,  and  phyllosilicates,  has  a  low  pressure  phase  characterized  by  po= 
2.615  Mg/m^-Kso  =  30.9  GPa,  K'  =  8.15,  and  y  =  0.73  (v/vo)'  «,  and  a  high- 
pressure  phase  with  po  =  3.659  Mgkm^,  Kjo  =  265  GPa,  K'  =  3.14,  and  y  =  1.2 
(VA^o)*  The  release  adiabats  calculated  from  these  equations  of  state  generally 
agree  with  the  data,  but  do  not  take  into  account  large-volume  phase  changes  that 
occur  at  low  pressure. 

Use  of  these  results  allows  us  to  estimate  the  fraction  f  of  shock  internal 
energy  irreversibly  deposited  in  these  rocks.  This  energy  in  no  longer  available  to 
drive  the  shock  wave,  so  the  shock  wave  decays  more  rapidly  than  expected  from 
geometric  effects  alone.  We  find  that  f  is  a  factor  of  ~2  smaller  in  limestones  than 
in  sandstone,  so  that  decoupling  of  the  source  from  far-field  seismic  waves  is  more 
efficient  in  sandstones.  This  analysis  ignored  the  effects  of  material  strength  at  low 
shock  pressures,  although  such  effects  are  important.  The  value  of  f  may  be 
reduced  considerably  in  the  presence  of  an  elastic  precursor  wave  and  of 
intermediate  waves.  Experimental  evidence  indicates  that  shock  wave  decay  still 
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occurs  even  in  the  supposedly  purely  elastic  shock  regime,  showing  that  effects 
such  as  fracture  formation  are  very  impoiiant  to  lire  overall  picture  of  explosive 
energy  coupling. 
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